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Abstract 

The analysis of nonlinear interaction of transversal electromagnetic field 
with collisionless plasma is carried out. Formulas for calculation electric current 
in collisionless plasma with arbitrary degree of degeneration of electronic gas 
are deduced. It has appeared, that the nonlinearity account leads to occurrence 
of the longitudinal electric current directed along a wave vector. This second 
current is orthogonal to the known transversal current, received at the classical 
linear analysis. 
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Introduction 


Dielectric permeability in quantum plasma was studied by many authors 
[i] - ra. Dielectric permeability is one of the major plasma characte¬ 
ristics. 

This quantity is necessary for the description of skin-effect [12], for the 
analysis surface plasmons [13], for descriptions of process of propagation 
and attenuation of the transversal plasma oscillations [8], for studying 
of the mechanism of penetration electromagnetic waves in plasma [7], 
and for the analysis of other problems in the plasma physics PI-[HI. 
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Let us notice, that for the first time in work [L] the formula for calcu¬ 
lation of longitudinal dielectric permeability into quantum plasma has 
been deduced. Then the same formula has been deduced and in work 


In the present work formulas for calculation electric current into 
collisionless plasma at any temperature (at any degrees of degeneration 
of the electronic gas) are deduced. 

It has appeared, that electric current expression consists of two sum¬ 
mands. The first summand, linear on vector potential, is known classical 
expression of electric current. This electric current is directed along 
vector potential of electromagnetic held. The second summand represents 
itself electric current, which is proportional to the square vector potential 
of electromagnetic held. The second current is perpendicular to the first 
and it is directed along the wave vector. Occurrence of the second current 
comes to light the spent account nonlinear character interactions of 
electromagnetic held with plasma. 

Let us underline, that the nonlinear phenomena in plasma are studied 
more half century (see, for example, [2D] - [22] )■ However, the electric 
current directed along a wave vector, was not is revealed earlier. 


1. Vlasov kinetic equation and its solution 


Let us consider Vlasov equation describing behaviour of collisionless 
plasmas 



(i.i) 


Vector potential we take as orthogonal to direction of a wave vector 



( 1 . 2 ) 


i(kr— ut) 


A(r, t) = A 0 e' 
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Scalar potential we will consider equal to zero. Electric and magnetic 
fields are connected with vector potential by equalities 

Id A 


E = 


c dt ’ 


H = rot A. 


(1.3) 


The wave vector we direct along axis x: k = fc(l,0,0), and vector 
potential of elec tromagnetic field we direct along axis y: 


A = Ay(x,t)( 0,1,0), Ay(x,t) ~ e^ kx ~ ut \ 


Then 


ic „ 

H = 

ck 

— Ey, 


uj y 

ck 

u 

[v,H] 

= 

v y, 


Let us operate with method of consecutive approximations. Consi¬ 
dering, that the member with an electromagnetic field has an order, on 
unit smaller other members, let us rewrite the equation (1.1) in the form 

q f( k ) d fW 

-0 v x - L 


dt 


dx 


+eiJ „|^Cq_h£U^^I=0. k = 1,2. (1.4) 


dp y 


U) 


lu dpj 


Here in zero approximation is the absolute Fermi Dirac distri¬ 
bution, 


/ (0) = /o = 


/ mv 


1 + CXP l2 ksf-^f)] 


k B is the Boltzmann constant, T is the plasmas temperature, p is the 
chemical potential of plasmas. 

It is easy to see, that 

.2 „ £ 


mv 


p 


=- a = P 2 — a, 

2 k B T k B T £ r 
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where a is the dimensionless (normalized) chemical potential, P = 
v/vt = p/pr is the dimensionless electron velosity (or momentum), 


£ = 


mv 


£r 


in r 


T 


2 7 2 
Therefore in zero approximation 


Vt 


'2k B T 


m 


a 


p 


k B T 


fo(P) = 


1 + e 


P 2 -c 


We notice that 


because 


V - H i ~ir = 0 ' 


dfW 

dp 


rv V. 


We search solution in first approximation in the form 

f m = f {0) + h(x,t,P x ). 

In this approximation the equation (1.4) becomes simpler 


<9/i dh 
+ v. 


= —eE, 


6>/(°) 


dt ' x ax 
From the equation (1.5) it is found 

2 ieE y P y g(P ) 


°Py 


fi = 


Pt u ~ kv T P x 


where 


x P 2 -a 


g(P) = 


(1.5) 


( 1 . 6 ) 


(l + e p2 -*) 2 ' 

In the second approximation for function /^ we search in the form 
f( 2 ) = j(i) + f 2 ( Xj t,v x ), where / 2 - E%(x,t). 

Let us substitute /^ 2) in (1.4). We receive the equation 

a/ 2 df 2 dfi e r TTl d/i 

+ W# + + - fv, HI-// = 0. 


at 


ax 


am 


c 


ap 
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From this equation it is found 

d 


f2 = 


e 2 El 


P 2 t {uj - kv T P x ) 


y 


dp y \ p y g(p) / | kv Tr , d 

u — kv T P x UJ 

d 'p v g(P)\ 


g(P ) 


dP x \uj — kv T P x 


kVT p OPy 


UJ 


uj — kvrP x 1 ^ ^ 

Distribution function in square-law approximation on the held it is 
constructed 

/ = / (2) = / (0) + fi + A, (1-8) 

where /i, are given accordingly by formulas (1.6) and (1.7). 

2. Density of electric current 


Let us calculate current density 

j = e [ vf 


2 d 3 p 


J pTTfiV (2 ' X) 

By means of (1.8) it is visible, that the vector of current density has 
two nonzero components 

j ( Jxijyi 0) • 

Here j y is the density of known transversal current, calculated as 

( 2 . 2 ) 


/ , 2 d 3 p f , 2 d 3 p 

h = e v v -(= e / v vfy 


(2irh)3 J yJ 1 (2tt h) 3 

This current is directed along electric held, its density is deduced by 
means of linear approximation of distribution function. 

Square-law on quantity of an electromagnetic held composed the 
contribution to density of a current does not bring. Density of transversal 
current it is calculated under the formula 


h ~ 


- 2 T EJx,t) l' P P P)d3P 


27r 3 m 


uj — kv'jP x 


(2.3) 
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Here kx is the thermal wave number, 


Pt ijivt 

T ~ h ~ h 

Let us calculate the longitudinal current. For density of longitudinal 
current according to definition it is had 

f 2 d 3 p f 2 d 3 p 

Jx = e I v x f /rt 3 = e / v x f r 


( 27 t ny 

2evxPx 


( 27 t ny 


P;j2d S P- 


(27tH) 3 j 

Having taken advantage (1.7), from here we receive, that 

dP, \ Py9 ^) kv T 2 d ( g{P ) 


_ 3 2 2mvj. 
h - V (27 Thy 


uj — kvrP x 

d 'Pyg(P) 


kv T „ dP< ^ y 


■Pr 


uj uj — kv T P 
Equality (2.4) can be simplified 

1 d 


uj y dP x \uj — kvrP x 


P r (PP 


uj — kvxP x 


(2.4) 


O O 2TTblJrp 
J X = PEl T 


(2ti ny 


ujdP y ~ y 


( Py9(P)) + 


kvT p 2 d 


9(P ) 


P r (PP 
uj — kvrP x 


uj ~ y dP x \uj - kv T P x 1 1 ■' 7 - n 

The first integral from (2.5) is equal to zero. Really, we will consider 
internal integral on P y 


OO 


/ w (Py9(F))dPy = Py9(p) 


— OO 


Py = + OO 


= 0. 


Py = — 00 


Now in the second integral from (2.5) we will calculate internal integral 


on P r 


OO 


d 


9(P) 


dP x \uj — kvrP x ' w — kv T P, 


P x dP x 


— OO 
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g(P)Px 


u - kv T P x )) 2 


P x =+00 00 


3(P) 


-d 


P =—oo _ 


cu — kvxP x — kvrPj 


P, 


OO 


oo 


—a; 


g(P)dP x 


(w - kv T P x ) 3 ' 

—OO 

Thus, equality (2.5) becomes simpler 


j* = —e E. 


5l7l 2^ mV T^ f 9(P)Pyd3P 


1:1 (27Th) 3 J (u - kv T P x ) 3 


e 3 El 


g(p)pyp 


PE 3 q f g(P)P 3 d 3 P 


\'K 3 kimv\q 2 J (P x — kl/q) 3 \.ix 3 kimv\ J (qP x — k}) 3 
Here 


ft = 


Cd 


kxVT 


k 

q = V T 


( 2 . 6 ) 


Double internal integral from (2.6) in plane (P y , P z ) it is calculated 
in polar coordinates (P y = pcostp , P z = psin(p): 


oo oo 27 t oo 2 

g(P)PydPydP z = J j -^^W~^T d( pdp = 

—oo —oo o o 

= ^ln(l + e“m )j P 3 = p + Pl 
Hence, the density of longitudinal current is equal to 

n o 00 

. _ e 3 E]q r ln(l + e a ~ Px )dP x __ 

4:7T 2 hmv^q 2 J (P c —ft/g ) 3 

— OO 



or 


oo 

e 3 E 2 q r ln(l + e a ~ p *)dP x 

47 i 2 hmv^ J ( qP x — ft) 3 

—OO 


e 3 E 2 


A^hmv^q 1 


J \n(l+e a - {p * +n/q)2 


dP x 

p3 ■ 


(2.7) 


( 2 . 8 ) 


—OO 



Let us present the formula (2.8) in an invariant form. For this purpose 
we will introduce transversal electric held 


„ „ k(Ek) „ q(Eq) 

J^tr I--TV- I-- 9 -• 

Ar q 

Now equality (2.8) we will present in coordinate-free form 


OO 


•long_ 

•3 


e 3 Ep.q f ln(l + e a T ~)dr 




(i qr — O) 3 


— OO 


The integral (2.7) is calculated with use of known rule of Landau as 
Cauchy type integral 


47 i 2 hmv^ 


_ CT J_( in(i + e «-yy 


+ 


T=fl/q 


+V.p. 


ln(l + e a t2 )cIt 
( qr — O) 3 


— OO 


Symbol V.p. means principal value of integral. The sign choice x — ie 
in the previous formula means attenuation in time of potential of an 
electromagnetic held 


Aoe z(kr^-ie)t) = AQe -et e i(kr-ut) > Q. 


We will present this formula in the following form 

j long = WtrE? r q./(0, q), 

where J(Q,q ) is the dimensionless part of electric current density 

ln(l + e a ~ T ~)dr 


J(a«) = -W|ln(l + e“-p 


1 rt 


OO 


+ V.p. 


r=tt/q 


—OO 


(qr — f2) 3 
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e 3 k B T 


Cltr = 


47 x 2 hmv\ 4:7T 2 hm 2 


On fig. 1-3 we will present graphics of behaviour of the real part of 
dimensionless quantity of density of electric current 


Re( J(0, q)) = V.p. 


ln(l + e a T ~)dr 
(<qr — O) 3 


— OO 


On fig. 4,5 we will present graphics of behaviour of the imaginary part 
of dimensionless quantity of density of electric current, 


Im (J (0, q)) = 


71 


ln(l + e a_r2 )j " 


7T 


2 q 3 L 

(1 - 2r 2 )e a ~ r2 + e^ a ~ T ^ 
q 3 ( 1 + e a ~ r2 ) 2 

1 + (1 - 2 T 2 )e r2 ~ a 


r=fl/q 


r=fl/q 


71 


r=Q/q 


q 3 ( 1 + e r2 - a ) 2 

In case of small values of wave number from (2.6) it is received 

o o 00 

e 3 E 2 (x,t) 


27 Phmv^Pl 3 


q ■ / ln(l + e a T )dr. 


It is possible to present this equality in the form 


J x = - 2e>E ^ t ^ T • Q ■ f ln(l + e a ~ r2 )dr. 


7i 2 h{hw) 3 

0 

Let us find numerical density (concentration) of plasma in the equi¬ 
librium condition 


N = J fo(P) 


2 d 3 p 87 rp: 


OO 0 

3 r e a-P 2 p 2 


P 2 dP 


(2nh) 3 (27ih) 3 J 1 + e a p2 

0 
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oo 

= ft / ln(1+ ea - p2)dR 

o 

Integral from expression for current density we will express through 
numerical concentration of plasma in equilibrium condition. It is as a 
result received, that 

eEy(x , t) ujp eEy(x , t ) uj 2 

^ x 47r (hco) UJ 2 ^ 47 IIJIVtU cu 2 

Here u p is the plasma (Langmuir) frequency, 


UJ r 


' Aixe 2 N 


m 


In coordinate-free form last equality rewritten as follows 



e K ", 

47r(hcc) c 2q 


3. Conclusions 

In the present work the solution of Vlasov equation is used for collisionless 
plasmas. For the solution it is used the method of consecutive approximations. 

As small parametre the quantity of the vector potential of electromagnetic 
held (or to it proportional quantity of intensity of electric held) is considered. 

At use of approximation of the second order it appears, that the 
electromagnetic held generates an electric current directed along the 
wave vector, and proportional to the size square of electric held. 

Thus an electric current directed along electric held, the same, as in 
the linear analysis. 
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Fig. 1. Real part of longitudinal electric current density, 0 = 1. Curves 1,2,3 
correspond to values of dimensionless chemical potential a = —1, —0.3, +0.1. 



q 


Fig. 2. Imaginary part of longitudinal electric current density, 0 = 1. Curves 1,2,3 
correspond to values of dimensionless chemical potential a = —1, —0.3, +0.1. 
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Fig. 3. Real part of longitudinal electric current density, a = 0. Curves 1,2,3 
correspond to values of dimensionless wave number q = 0.5, 0.4, 0.3. 



Fig. 4. Imaginary part of longitudinal electric current density, a = 0. Curves 1,2,3 
correspond to values of dimensionless wave number q = 0.5, 0.4, 0.3. 
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Fig. 5. Real part of longitudinal electric current density, a = 0. Curves 1,2,3 
correspond to values of dimensionless frequency of electromagnetic field 

Q = 1,1.2,1.5. 
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Fig. 6. Imaginary part of longitudinal electric current density, a = 0. Curves 1,2,3 
correspond to values of dimensionless frequency of electromagnetic field 

O = 1.1.2.1.5. 











